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ABSTRACT 

Quasars (QSOs) at the highest known redshift {z ^ 6) are unique probes of the early growth of 
supermassive black holes (BHs). Until now, only the most luminous QSOs have been studied, often one 
object at a time. Here we present the most extensive consistent analysis to date of 2: > 4 QSOs with 
observed NIR spectra, combining three new z 6 objects from our ongoing VLT-ISAAC program with 
nineteen 4 < z < 6.5 sources from the literature. The new sources extend the existing SDSS sample 
towards the faint end of the QSO luminosity function. Using a maximum likelihood fitting routine 
optimized for our spectral decomposition, we estimate the black hole mass (Mbh), the Eddington 
ratio (defined as Lboi/LEdd) and the Fell/Mgll line ratio, a proxy for the chemical abundance, to 
characterize both the central object and the broad line region gas. The QSOs in our sample host BHs 
with masses of ^ 10^ that are accreting close to the Eddington luminosity, consistent with earlier 
results. We find that the distribution of observed Eddington ratios is significantly different than that 
of a luminosity-matched comparison sample of SDSS QSOs at lower redshift (0.35 < z < 2.25): the 
average {log{Li)oi/ LEdd)) = —0.37 {Li^oi/LEdd ~0.43) with a scatter of 0.20 dex for the z > 4 sample 
and the {log{Lboi/ LEdd)) = —0.80 {Lijoi/LEdd ~0.16) with a scatter of 0.24 dex for the 0.35 < z < 2.25 
sample. This implies that, at a given luminosity, the Mbh at high-z is typically lower than the average 
Mbh of the lower-redshift population, i.e. the z > 4 sources are accreting significantly faster than the 
lower-redshift ones. We show that the derived Fell/Mgll ratios depend sensitively on the performed 
analysis: our self-consistent, homogeneous analysis significantly reduces the Fell/Mgll scatter found 
in previous studies. The measured Fell/Mgll line ratios show no sign of evolution with cosmic time 
in the redshift range 4 < z < 6.5. If the Fell/Mgll line ratio is used as a secondary proxy of the 
Fe/Mg abundance ratio, this implies that the QSOs in our sample have undergone a major episode of 
Fe enrichment in the few 100 Myr preceding the cosmic age at which they are observed. 
Subject headings: cosmology: observations - quasars: general, emission lines - galaxies: active, high- 
redshift, formation 



1. INTRODUCTION 



In the past 10 years, more than 50 QSOs a t z > 



5.7 have b e en discovered (iFan et al.' '200 d |2QQli 
2OOI [2006[ iJiang et al.l [20(^ 2009: Willo tt et all 



2003, 



2007, 



2OIOD thanks to the Sloan Digital Sky Survey (SDSS 



York et al]l2000f ) and other large multi- wavelength sur- 



veys such as the Canada F rance High-z Quasar Sur- 
vey fCFHOS. IWillott et al.l 12007). These high-redshift 
QSOs are among the most luminous sources known to 
date and are direct probes of the universe less than 1 
Gyr after the Big Bang. They are fundamental force 
studying the early growth of supermassive BHs, galaxy 
formation and interstellar mediu m chemical evolution 
(e.g. Kauffmann Haeneltl 120001 : IWvithe L"Qebl 120031 : 

iHopkins et al. 20051 ) . Strong emission features excited by 

the central engine and its immediate surroundings can be 
used to infer properties of the powering BH and of the 
circumnuclear gas. For example, measuring the emission 
line width one can obtain information about the bulk mo- 
tion of the broad-line region (BLR) gas, that can then 
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be used to estimate the BH mass (Mbh)- From the line 
flux ratios we can instead derive chemical abundances of 
the BLR gas that are fundamental to set constraints on 
the star formation history of the QSO host galaxy. 

Under the assumption that the dynamics of the BLR 
is dominated by the central gravitational field, it is pos- 
sible to estimate Mbh by using the velocity and the 
distance of the line-emitting gas from the central BH 
: Mbh oc R v'^. From reverberation mapping studies of 
local active galactic nuclei (AGNs) it has been found that 
R is related to the continuum luminosity L {R — L re- 
lation, e.g. iKaspi et al.|[2QQQf ): R oc L^-^. Measurements 
of the luminosity and of the gas velocity, via Doppler- 
broadened line-width, can thus be used to determine the 
Mbh from single epoch spectra. Various strong emission 
lines can be used as Mbh estimators: Ha, H/3, Mgll, 
C I V (for a review on the single-epoch spectrum method 
see iPetersonI |2Q1Q[ ). For objects in the local universe, 
the R — L relation and its intrinsic scatter have been in- 
vestigated in detail using the H/3 line and the relative 
continuum. These studies have shown that one can ob- 
tain accurate Mbh estimates via the H/3 line. For sources 
at high-z a complication is that the H/3 emission line is 
redshifted out of the visible window already at modest 
redshift s. Whether or not a particular UV line can be 
used to estimate Mbh depends on how well the respec- 
tive UV line wi dths in AGN spect ra are correlated. In 
the case of Mgll lShen et al.l (|2QQ8[ ). using a large collec- 
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tion of SDSS spectra, found that log ([FWHM(H/3)] / 
[FWHM(Mgll)]) = 0.0062, with a scatter of only 0.11 
dex, suggesting that the Mgll can be used as a proxy 
for the H/3 line, and thus for the BH mass. On the 
other hand, the accuracy achievable in the determina- 
tion of the BH mass from the C IV line remains some- 
what controversial, since C IV is a resonance line and 
absorptions are often detected in its blue wing due to 
outflows, which in turn affects the line width. With the 
single-epoch spectrum method it has been possible to 
estimate the Mbh for z - 6 QSOs (Barth et al. 20031; 
IWillott e t all 120031: iJiang et al.l 120071: iKurk et al.l I2007L 
|2009; W illott et al.l 12010^. These studies have shown 
that these high-z sources from SDSS host BHs with Mbh 
^ 10^ M0 and are accreting close to the Eddington limit 
{Li)oi/LEdd "^1)- However, the z ^ 6 QSOs from SDSS 
only account for the brightest end of the QSO luminosity 
function. Only few faint QSOs (^^^ > 20) have mea- 
sured Mbh to date: two f r om the SPSS deep Str ipe 
82 Jiang et al.l [20081 [2QQ1 iKurk et al.l [2QQ7L iSoOOl ). a 
deep imaging survey obtained by repeatedly scanning a 
stripe (260 deg^) along the celestial equator, and nine 
from CFHQS (Willott et al. 2010). These sources are 
powered by less massive BHs (down to Mbh ^ 10^ ^o) 
that also appear to accrete close to the Eddin gton limit. 
For QSOs at lower redshifts lShen et al.l (I2008D found that 
the most luminous QSOs (L^oi > 10^^ erg s~^) at red- 
shift 2 < 2; < 3 are characterized by an Eddington ratio 
of only Lijoi/LEdd ~0.25 with a dispersion of 0.23 dex. It 
seems therefore that high-redshift QSOs have fundamen- 
tally different properties than the lower-redshift ones. 

Regarding the chemical enrichment, photoionization 
models show that various emission line ratios are good 
estimators for the metallicit y of the BLR gas (e.g. 
iHamann et ani2002l : iNagao et" all 120061 ). It is possible to 
estimate the BLR gas metallicity and to set constraints 
on the BLR enrichment history through measurements of 
the relative abundances of nitrogen (N) respect to car- 
bon (C) and helium (He), since these three elements are 
formed by differ ent astrophysical processe s and on differ- 
ent timescales (jHamann fc Ferlaiidlll993f ): C is rapidly 
produced in the explosion of massive stars; N is a sec- 
ond generation element, i.e. slowly produced in stars 
from previously synthesized C and O; He is a primordial 
element and its abundance does not change significantly 
with cosmic time. Using these element ratios IJiang et al.l 
f2007) estimated the BLR metallicity for a sample of six 
luminous QSOs with 5.8 < z < 6.3, finding super solar 
metallicities with a typical value of ~ 4 Zq and no strong 
evolution up to z ~ 6. 

The abundance of Fe and a elements (e.g. O, Mg, Ne 
that are produced in the stars via a processes) is of par- 
ticular interest for understanding the chemical evolution 
of galaxies at high-z. Most of Fe in the solar neighbor- 
hoods has been produced via the explosion of type la 
supernovae (SNe la) while a elements such as Mg and O 
are mainly produced by core collapse supernovae (SNe of 
types II, lb and Ic). SNela are thought to originate from 
intermediate- mass stars in close binary systems, charac- 
terized by long life-times, while core collapse SNe come 
from more massive stars which explode very soon after 
the initial starburst. A time delay between a elements 
and the Fe enrichment is thus expected. This delay de- 



pends on the Initial Mass Function (IMF) and on the 
galactic star formation history, and can vary from 0.3 
Gyr, for massive elliptical g alaxies, to 1-3 Gyr, for Milky- 
Way type galaxies (Matteu cci fc Recchil l2003f ). In this 
picture the Fe/a ratio is expected to be a strong function 
of age in young systems. The observational proxy that is 
usually adopted to trace the Fe/Mg abundance ratio for 
the BLR gas is the Fell/Mgll line ratio. Unfortunately 
this line ratio is only a second order proxy, since it de- 
pends not only on the actual abundance of Fe but also 
significantly on the excitation conditions that determine 
how strongly Fe II lines are emitted (Baldwin et al. 2004). 
However, even though to date no calibration is available 
to convert the observed relative line strengths into actual 
abundance ratios, the study of the Fell/Mgll line ratio 
as a function of look-back time does in itself carry sig- 
nificant information about the BLR chemical enrichment 
history. 

For z > 5.7 QSOs the Fell and the Mgll lines 
are redshifted into the near-infrared (NIR). Numer- 
ous NIR-spectroscopy studies of QSO samples includ- 
ing z ^ 6 sources have been carried out in the 
past (e.g. Maiolino et al. 2001, 2003; Penteri cci et al.l 
20021: Iwamur o et al.l I2002L [20041: Barth et aD 120031: 
Dietrich et al. 20031: iF reudling et al.L2003: Willott et al.l 
2003; Jiang et al. 2007: iKurk et al.l 120071 . I2009f ). Their 
results show an increase in the scatter of the measured 
Fell/Mgll line ratios as a function of the redshift. One 
possible way to interpret the increase in the scatter is 
that some objects are observed such a short time after 
the initi al starburst that the B LR is not yet fully enriched 
with Fe (|Iwamuro et al.|[2004f ). Nonetheless, several pos- 
sible systematics related to the adopted fitting procedure 
could potentially be the cause for the observed discrep- 
ancies between different authors: the iron template em- 
ployed, the wavelength limits over which the template 
is integrated, the wavelength range over which the tem- 
plate is actually fitted and the wavelength coverage and 
S/N of the spectra used. 

The goal of this paper is to study BH masses, Edding- 
ton ratios and to characterize the Fe ll and Mg ll emission 
in a large sample of high-z QSOs, by fitting their respec- 
tive NIR spectra in a coherent and homogeneous way. 
With the aim of extending the existing SDSS sample of 
z ~ 6 QSOs towards the faint end of the QSO luminosity 
function, we observed three additional z ^ 6 sources (dis- 
covery papers: Jiang et al. 2008, De Rosa et al. in prep) 
with ISAAC (Infrared Spectrometer And Array Cam- 
era) mounted on the Very Large Telescope (VLT). The 
three targets have 19.6 < < 20.7 and have been 

selected from the SDSS and the SDSS Stripe 82. The 
observed NIR spectra include the Mgll and Fell emis- 
sion lines. We have also collected literature spectra of 
z > 4 QSOs covering the restframe wavelengths between 
2700 A< A < 3200 A, characterized by the presence of 
the Mgll line doublet (A = 2796, 2803 A) and of the Fell 
line forest. 

The paper is structured as follows. In § 2 we describe 
the sample, the new observations and the data reduction 
we have performed. In § 3 the spectral decomposition 
and the fitting procedure adopted are explained together 
with their results. In § 4 we estimate the Mbh and the 
Fell/Mgll line ratio and we discuss their implications. 
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Finally, we give a brief summary in § 5. We assume 
the following ACDM cosmology throughout the paper: 
Ho 70 km s-^ M pc~\ Qm 0.3, and Qa 0.7 
(jSpergel et al.ll2007[ ). 

2. SAMPLE, OBSERVATIONS AND DATA REDUCTION 

Our sample is composed of 22 targets: 3 new spectra 
observed with VLT-ISAAC (see Tab. H]) and 19 sources 
from the literature (see Tab. [2]), kindly provided by the 
respective authors. Ten of the literature so urces have 
redshift 4.50 < z < 5.70 (llwamuro et al.ll2002[ ). while the 
remaining 9 have 5.70 < z < 6.43 and z-band magnitudes 
^AB ^ 20-9 (magnitudes are taken from the discovery 
papers). 

2.1. New data 

We have observed 3 SDSS QSOs with magnitudes 
19.6 < z'^j^ < 20.7 and redshifts 6.05<z<6.08 (see 
Tab. [1]). The two faintest ones have been selected 
from the SDSS Stripe 82 and extend the existing sam- 
ple towards the faint end of the QSO luminosity func- 
tion. The observations were carried out with ISAAC 
on Antu (VLT-UTl) in low resolution mode (LR), using 
the 1024x1024 Hawaii Rockwell array of the Short Wave- 
length arm. For each QSO the Mgll line and the Fell 
complex were observed: given the redshift of the sources, 
these features fall in the K band. The selected slit had 
a width of 1" and combined with the order selection fil- 
ter it gives a spectral resolution A/ A A ~450. Table [1] 
summarizes the exposure time for each object. For each 
observation block (OB), sixteen frames of 148 seconds 
were taken following an ABBA dithering pattern (with 
large offsets among the dithered positions: from 20" to 
30"). Further, small random offsets within a box of 4" 
to 18" were applied at each dithered position in order to 
avoid pixel-related artifacts (jittering). Given the faint- 
ness of the sources, the observation setup was chosen so 
that a bright star was always in the slit, in order to allow 
a correct centering of the target. 

The QSO SDSS J2054- Q0Q5 has bee n discovered in the 
SDSS deep Stripe 82 bv iJiang et a T (2008). Our new 
ISAAC spectrum confirms the weak-line nature of this 
source (see also the the optical discovery spectrum by 
iJiang et al.l 120081 ). Given the intrinsic weakness of the 
Mgll emission, we do not include this QSO in the fol- 
lowing analysis. 

2.1.1. Data reduction 

The ESO ISAAC pipeline produces wavelength cali- 
brated co-added 2-D spectra from the individual frames 
that were acquired during each OB. Subsequent reduc- 
tion was carried out within IRAF. One-dimensional spec- 
tra were extracted using the apall task. The tracing of 
the 1-D spectra was performed first on the bright stars 
in the slit. The resulting tracing functions were then 
used for the extraction of the QSO spectra. Individual 
ID spectra were corrected for the telluric absorptions us- 
ing the telluric task. Usually the telluric correction is 
performed by dividing the observed spectrum by the one 
of a telluric standard star observed shortly after the sci- 
ence target. This ratio is subsequently multiplied by the 
model atmosphere corresponding to the spectral type of 
the telluric standard and scaled to its observed K mag- 
nitude, in order to recover the correct slope of the QSO 



spectrum. This operation allows us also to flux- calibrate 
the QSO spectrum. Instead of using the observed tel- 
luric standard stars, we employed the ESO sky absorp- 
tion spectrum measured on the Paranal site at a nominal 
air mass of 1. This choice was driven by two reasons: i) 
the spectral regions in the QSO spectrum where telluric 
absorptions are most severe are characterized by a very 
low signal-to- noise, i.e. insufficient to properly correct 
for the actual detailed shape of the night sky features; ii) 
the observed telluric standard stars were often character- 
ized by a spectral and luminosity class for which accurate 
model atmosphere could not be computed; the derived 
QSO continuum slope would have hence been distorted 
by the stellar spectral shape. We thus decided to focus 
on the spectral regions with higher signal, where an ac- 
curate correction was possible: we assumed the template 
sky absorption spectrum, scaled to the airmass of our 
QSO spectra and ignored the stability of the sky trans- 
parency. This way, we could preserve the intrinsic shape 
of the QSO continua. The telluric task corrects for the 
difference in airmass between the science and calibration 
spectra via the Beer-Lambert-Bouguer law. The telluric 
absorptions are in any case well removed from the spec- 
tra of bright QSOs, while significant residuals are left 
in the spectra of faint QSOs at lower S/N. The relative 
flux calibration was obtained with the sensf unction, 
standard and calib tasks. The instrument sensitivity 
function was obtained from the observed telluric stars 
of luminosity class V (giants) and spectral class B. For 
these stars it is possible to compute reliable model atmo- 
spheres since effective temperature and surface gravity 
are well estimated. The model atmospheres were com- 
puted by interpolating the N IR spectral energy distribu- 
tions available from P icklesI ([l998) in temperature and 
surface-gravity in order to match the observed spectral 
type and magnitudes. After the relative flux calibration, 
the individual ID QSO spectra were averaged to form a 
single spectrum. The absolute flux calibration was per- 
formed scaling the observed spectra to match the QSOs 
K-band magnitudes. Since the K-band magnitude was 
not measured for the QSOs in Tab. [TJ we derived it from 
the SDSS QSO template scaled to the observed J and 
H-band magnitude. The whole reduction procedure was 
extensively tested on some of the observed telluric stan- 
dard stars used as reference targets. The recovered spec- 
tra match the theoretical model atmospheres typically 
within 10%, even in the spectral regions more affected 
by telluric absorptions. The reduced spectra are shown 
in Fig. [331 

2.2. Literature Data 

The QSOs taken from the literature are summarized 
in Tab. [2l We have collected a total of 22 spectra (19 
different sources) which cover all the features of inter- 
est at sufficiently high S/N to perform the spectral de- 
composition (see Sec. 14.21 and Sec. 14. 6p . The literature 
sample is composed of 10 sources with z < 5.70 and 
9 sources with 5.70 < z < 6.43. The 10 QSOs with 
z < 5.70 are selected am ongst the 13 sources published 
by llwamuro et al.l (|2002f ). These QSOs have redshifts 
4.4 < z < 5.3 and were observed in the J and H bands 
with an OH-airglow suppressor spectrograph (OH-S), 
mounted on the Subaru Telescope. The spectral resolu- 
tion (A/AA) is equal to 210 (J band) and 420 (H band). 
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TABLE 1 

Summary of the newly observed sources. (1) QSO name; (2), (3) QSO right 

ASCENSION AND DECLINATION IN J2000.0 (fROM DISCOVERY PAPERS, COL [8]); (4) OPTICAL 

REDSHIFT (from DISCOVERY PAPERS); (5) z\q (fROM DISCOVERY PAPERS); (6) K BAND 
NOMINAL MAGNITUDE, OBTAINED FROM THE REDSHIFTED SDSS QSO TEMPLATE SCALED TO 
MATCH THE OBSERVED J AND H-BAND MAGNITUDES (DISCOVERY PAPERS); (7) TOTAL ON 
SOURCE EXPOSURE TIME; (8) DIS COVERY PAPER : DRiP , De ROSA ET AL. (iN PREP); J08, 

IJlANG ET AlI (l2008l) . 



QSO name 


RA 


DEC 


z 




K 


texp 


Reference 


(1) 




(2) 


(3) 


(4) 


If 


(6) 


(7) 


(8) 


SDSS J0353 + 


0104 


03 53 49.72 


+01 04 04.4 


6.05 


20.5 


17.81 


4.60 h 


J08 


SDSS J2054 - 


0005 


20 54 06.49 


-00 05 41.8 


6.06 


20.7 


17.87 


5.92 h 


J08 


SDSS J0842 + 


1218 


08 42 29 


+12 18 50.5 


6.08 


19.6 


17.67 


4.60 h 


DRip 



We also selected 2 of the 4 SDSS QSOsatz ~ 6 presented 
by llwamuro et al.l ()2QQ4f ): H and K observations were 
carried out with the Cooled Infrared Spectrograph and 
Camera (CISCO) mounted on the Subaru telescope, with 
a spectral resolution of 210 and 330 respectively. We used 
the K-band s pectrum of t h e z = 6.43 QSO J1148+5251 
published by IB art h et al.l (|2QQ3l ). Data were obtained 
with the NIRSPEC spectrograph on Keckll (spectral res- 
olution of 1500). We selected 4 of the 5 sources pre- 
sented by iKurk et all (I2007D at redshifts z > 5.8 . The 
K-band obseryations were carried out with VLT-ISAAC 
with a spectral resolution of 450. We also i ncluded the 
K-ban d spectrum of J0303-0019 presented by lKurk et al.l 
(200§). This faint z - 6 QSO was selected in SDSS 
Stripe 82 and its spectrum was taken with VLT-ISAAC. 
Finally, we added 4 SDSS QSOs at z > 5.8 published by 
iJiang et al] (|2P07), who obseryed them with GNIRS on 
Gemini South with a spectral resolution of 800. 

3. DATA ANALYSIS 
We focus on the spectral region with restframe waye- 
lengths 2000 A < Xrest < 3500 A. This region is charac- 
terized by the presence of the Mg ll emission line, the un- 
derlying non-stellar continuum, the Balmer pseudo con- 
tinuum and the Fell emission line forest. The last three 
emission features, that are oyerlapped in the spectral 
range of interest, are described in the following sections. 

3.1. Power-law 

The dominant component of a QSO spectrum is the 
non-stellar continuum, modeled as a power-law: 

Typically, in our data, the determination of the slope 
coefficient a depends on the adopted fitting procedure 
and on the obseryed spectral range. In case of a wide 
wayelength coyerage it is possible to choose fitting win- 
dows free of contributions by other emission compo- 
nents. For spectra with restricted wayelength coyerage, 
the power-law continuum has to be fitted simultaneously 
with the other components, resulting in a local estimate 
of the slope that might not be fully representative^^ 
the ov erall continuum shape of the QSO. iDecarli et al.l 
(|2010[ ) analyzed a sample of 96 QSOs at 2: < 3 by fitting 
the power-law continuum in 8 different windows free of 
strong features, and obtained a mean valu e for th e slope 
of —1.3 with a 1-cr dispersion of 1.6. .Shen et al.l (j2010.) 



estimated the local slope for a sample of ~ 100.000 SDSS 
QSOs at z < 4.95, fitting the the power-law plus an 
iron template to the wavelength range around four broad 
emission lines (Ha, II/3, Mgll and Civ). In particular, 
from the analysis of the wavelength range adjacent the 
Mgll emission line, they obtained a mean value for the 
local slope of —1.3 with a 1-a dispersion of 0.8. From 
our spectral decomposition we find consistent values (see 
Sec.lMl). 

3.2. Balmer pseudo- continuum 

We model t he Ba lmer pseudo-continuum following 
iDietrich et all (j2003f ). We assume partially optically 
thick gas clouds with uniform temperature Tg = 15000 
K. For wavelengths below the Balmer edge {Xbe = 3646 
A), the Balmer spectrum can be parametrized as: 

Fx = F„Bx{T,){l-e-^^^^^^'), A < \be (2) 

where Bx(Te) is the Planck function at the electron tem- 
perature Tg, tbe is the optical d epth at the Balme r edge 
(we assume tbe = 1 following iKurk et al.l 120071 ). and 
Fnorm is the normalized flux density at the Balmer edge 
(Grandi 1982). The normalization should be determined 
at Xrest — 3675 A, where no Fell emission is present. 
Since this wavelength is either not covered or has a very 
low S/N in our sample, we fix the normalization to a 
fraction of the continuum strength extrapolated at 3675 
A: Fnorm = fs ' Fpower-iawi^Q^^ A). To define the rela- 
tive strength of the two components and monitor the ef- 
fects on the Fe ll estimate, we have run various tests with 
fs = 0.1,0.3,0.5,0.8,1. Since differences in the Fell es- 
timates resulting from this test were less than measured 
errors (with /b only partially affecting the power-law 
normali zation), we have fixed /b = 0.3 based on the re- 
sults bv IDietrich et al.l (|2003D . 

3.3. Fell template 

The Fe ll ion emits a forest of lines, many of which are 
blended. We fit the Fell forest us in^ a modified version 
of the emission line template by 'Vestergaard fc WilkesI 
(|2001f ). This template is based on the high reso- 
lution spectrum of the narrow line Seyfert 1 galaxy 
PG0050+124 (z=0.061), observed with the Hubble Space 
Telescope. Vestergaard & Wilkes obtained the emission 
line template by ffist fitting and subtracting the power- 
law continuum and all the absorption/emission features 
from all the elements but Fe. Afterward an Fe ill model 
was subtracted from the residual to obtain a pure Fell 
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TABLE 2 

Literature sample sorted by redshift. (1) QSO name; (2), (3) QSO right 

ASCENSION AND DECLINATION IN J2000.0; (4) OPTICAL REDSHIFT FROM DISCOVERY 

PAPER; (5) NIR SPECTRUM REFERENCE PAPER FOR THE NIR SPECTRA: A 
[H^MURO ET AL. (2002), B IIWAMURO ET A L I (120041) . C IBarth et al.i (i2003i) . D 
IKURK ET A L. (2007), E IKURK ET AL.I (120091) F IJlANG ET AL.I (120071) . 



(^dU name 




RA 




z 


Reference for NIR spectra 






(2) 




(4) 




BR 1033 - 0327 


10 


33 


51.47 


-03 27 45.5 


4.51 


A 


BR 0019 - 1522 


00 


19 


35.90 


-15 22 16.0 


4.53 


A 


BR 2237 - 0607 


22 


37 


17.44 


-06 07 59.7 


4.56 


A 


SDSS J0310 - 0014 


03 


10 


36.97 


-00 14 57.0 


4.63 


A 


SDSS J1021 - 0309 


10 


21 


19.16 


-03 09 37.2 


4.70 


A 


SDSS J0210 - 0018 


02 


10 


43.17 


-00 18 18.4 


4.77 


A 


SDSS J0211 - 0009 


02 


11 


02.72 


-00 09 10.3 


4.90 


A 


PC 1247 + 3406 


12 


47 


17.79 


+34 06 12.7 


4.90 


A 


SDSS J0338 + 0021 


03 


38 


29.31 


+00 21 56.3 


5.00 


A 


SDSS J1204 - 0021 


12 


04 


41.73 


-00 21 49.6 


5.03 


A 


SDSS J0005 - 0006 


00 


05 


52.30 


-00 06 56.0 


5.85 


D 


SDSS J1411 + 1217 


14 


11 


11.30 


+12 17 37.0 


5.93 


D, F 


SDSS J1306 + 0356 


13 


06 


08.20 


+03 56 26.0 


5.99 


D, F 


SDSS J1630 + 4012 


16 


30 
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template. As iKurk et al.l (j2QQ7[ ) pointed out, no Fell 
emission is left in the 2770-2820 A range, due to the 
Mg II line sub traction. We modified the template follow- 
ing K urk et a l. (2007) by adding a constant fiux density 
between 2770 and 2820 A equal to the 20% of the mean 
fiux density of the template between 2930 and 2970 A. 
The justification for this operation co mes from the the- 
oretical Fell emission line strength by Si gut fc Pr adhanI 
(2003). Since the Fe ll and the Mg ll are not distributed 
in the same way within the BLR, we decided not to 
fix the Doppler broadening of the Fell template to the 
one measured for the Mgll emission line. We have 
instead run many tests broadening the Fell template 
by convolving it with Gaussian profiles with constant 
FWHM = 7.5, 15, 22.5 A, corresponding to FWHM 
- 930, 1860, 2800 km/s at 2400 A respectively. A Gaus- 
sian broadening with a constant in wavelength leads to 
slightly different velocities over the fitting range, but we 
have found no significant differences in the measured Fe ll 
normalization in function of the three broadening. Given 
the typical S/N of our spectra the Fe ll normalization is in 
fact mainly determined by the fit of the template broad 
bumps rather than individual Fell features. We finally 
chose to fix FWHM= 15 A for the Fell template. 

3.4. Fitting procedure 

All the spectra have been shifted to the rest-frame sys- 
tem of reference using optical redshifts. Even if the opti- 
cal redshift s are slightly different from the NIR ones, this 
will not affect our results since the Mg ll peak wavelength 
is a free parameter in our fitting procedure. Moreover, 
the Fell template is constituted by blended multiplets 
(broadened by the convolution with Gaussian profiles 
with FWHM = 15 A), for which the definition of peak 
wavelength is not straightforward. We tested a posteriori 
the difference between the optical redshifts (see Tab. [1] 
and Tab. [2j) and the ones measured from the Mgll peak 



wavelength (see Tab.|4]). The average difference is equal 
to Az=0.02 with a maximum of Az=0.08 for J1411 + 1217 
(j Jiang et al.ll2007L online Appendix Fig. A. 12) : in this 
case the Mgll line profile is severely affected by the at- 
mospheric absorptions. Our spectra have different wave- 
length coverage because of the different redshifts of the 
sources and of the various instruments used to collect the 
data. Similarly the sky contamination varies across the 
sample with redshift. For these reasons it is impossible 
to define fixed fitting windows for the entire sample. We 
focus on the rest-frame region within 2000 and 3500 A, 
choosing the fitting windows in a way as homogeneous as 
possible, as a function of the spectral coverage and of the 
sky contamination. The fit is performed in two steps. A 
first set of spectral components is given by the sum of the 
power-law continuum, the Balmer pseudo-continuum and 
the Fell template. The free parameters are the power- 
law slope {a) and its normalization {/3 = log{Fo)) and the 
Fell template normalization (7). Since the three compo- 
nents are overlapped in this wavelength range, the fit is 
performed with a minimization on a suitable grid in 
the parameter space. In this way, we minimize the possi- 
bility that our solution represents only a local minimum 
of the domain. The errors on the Fell normaliza- 
tion are computed by marginalizing the probability dis- 
tributions in 3-D parameter space, selecting all the cases 
for which ~ Xmin < 1 (IcJ" confidence level). To en- 
sure a reliable estimate of these errors we need enough 
triplets satisfying the condition. After many tests we 
decided to sample the parameter space with 16 million 
points: 100 (a) x 100 {(3) x 1600 (7). This way, for 
each fit, we sample the 7 parameter space with an av- 
erage bin of 0.001 x 10~^^ erg s~^ cm~^ and we obtain 
~ 100 triplets satisfying the with a minimum of 15 
(for the QSO BR 1033 - 0327). This cond ition is not 
satisfied for J0303 - 0019 (|Kurk et all 120091 . online Ap- 
pendix Fig. A. 17): in this case the Fe ll major features are 
severely affected by the telluric absorptions, resulting in 
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a weaker constraint on the Fell template normalization 
which is in any case consistent with (see Tab. [3]). In 
Fig. 13. 41 we show two examples of the 3D x^-cube projec- 
tions and of the relative probability distribution for the 
Fe II template normalization (similar plots for all sources 
discussed here are shown in the online Appendix). The 
X^-maps are overall regular (non-patchy), implying the 
absence of secondary local minina. The degeneracy be- 
tween power-law slope and intercept is evident from the 
bottom-right plots in the two panels. 

The fitted components are then subtracted and we pro- 
ceed to fit the Mgll emission line. The Mgll line fit 
is performed with a least-squares procedure. Since the 
Mg II doublet is not resolved in the majority of our spec- 
tra, we model the emission line as a simple Gaussian 
(three free parameters: central wavelength, width and 
normalization). If there was a significant narrow second 
component that we do not resolve, this would lead to a 
slight underestimate of the black- hole mass. 

Examples of the spectral decomposition are shown in 
Fig. 13.41 The results of the fit and relative x^ maps for 
the literature sample are shown in the online Appendix. 
The fitted parameters are listed in Tab. [3l Even if we 
do not overcome the degeneracy between the power-law 
slope and intercept, we can compare the distribution of 
our local slope estimates with those in the literature. 
We obtain a mean value a = —1.5 ±1.2 which is in 
agreement within the uncert ainties with both the local 
slope estimate by Shen et^l] ()2008l ) and with the global 
one bv iDecarh et al.. (^20 IC 



4. RESULTS 

4.1. Black hole mass estimate 

We estimate the Mbh using scaling relations, cali- 
brated on local AGNs, that are based on broad emission 
line widths and continuum luminosities. Under the as- 
sumption that the dynamics of the broad line region is 
dominated by the gravity of the central BH, the virial 
theorem states: 



Mbh = G ^ Rblr v^lr 



(3) 



where Mbh is the black hole mass, Rblr is the charac- 
teristic radius of the BLR and vblr is the orbital veloc- 
ity of the clouds emitting at Rblr- The cloud velocity 
can be obtained from the width of the broad emission 
lines: vblr = / • FWHM, where / is a geometrical fac- 
tor that account s for the de-project i on of vblr, from th e 
line of sight fsee'McG ih et al.ll2"Q08l : iDecarli et aP 12008V 
and FWHM is the full width at half maximum of the 
line profile. Even if the BLR size cannot be directly 
measured by single epoch spectra, it can still be evalu- 
ated since it is strongly correlated wit h the con tinuum 
luminosity of the AGN (e.g. Kaspi et^ l2000[ ). It is 
then possible to estimate Mbh for high-redshift QSOs 
using a single spectrum covering the blended Mgll line 
doublet (A = 2795, 2803 A) and the redward continuum 
(A = 3000 A ). In particular for Rblr we use the relation 
provided by iMcLure fc JarvisI ()2002f ) : 



RBLRjMg II) 

10 light — days 



(2.52 ±0.03) 



ALa(3000 A) 



_ 10^4 grg s-i 



0.47 



(4) 



and for the geometrical factor the value provided by 
iDecarh et anl2008f ): /(Mgll) = 1.6, obtained assuming 
that the Mgll and H/3-emitting regions have a similar 
geometry. All these relations are based on low redshift 
objects {z < 0.3), and the underlying assumption here is 
tha t they are still va li d at h i^h-redshift. Using relation 
dH) iMcLure fc JarvisI (|2002r ) reproduced Mbh obtained 
with the full reverberation mapping method with an ac- 
curacy of 0.4 dex. This intrinsic scatter of the estimator 
dominates the measurement uncertainties. To compare 
our results with the ones published in previous works 
on z - 6 sources ( Barth e t all 120031 : iJiang et al.l 120071 : 
iKurk et al.ll2007L 12009). we also estimat e Mbh using; the 
relation obtained by iMcLure fc Dunlopl (|2C 



M 



BH 



Mr. 



3.2 



ALa(3000 A) 
10^7 W 



0.62 



FWHMiMg II) 



km s 



(5) 

from a sample of 17 low redshift AGNs {z < 0.7), with 
luminosities comparable to those of high-z QSOs {XLx > 
10^"^ erg s-^). 

From the estimated Mbh we compute the QSO Ed- 
dington ratios defined as the ratio between the measured 
bolometric luminosity Lijoi and the theoretical Eddington 
luminosity LEdd^ as computed from the Mbh- We obtain 
the observed monochromatic luminosity ALa(3000A) 
from the continuum component of the fit /a (3000 A) and 
from the luminosity distance Dl as: 



ALa(3000A) = A Fa (3000 A) 47rD^ 



(6) 



Since L (3000 A) is only a fraction of the total electro- 
magnetic luminosity comin g from t he QSO, we apply 
the bolometric correction bv lShen et al.l ((20081) to obtain 
Lboi- 

Lboi = 5.15 ALa(3000A) (7) 

The Eddington luminosity is defined as the maximum 
luminosity attainable, at which the radiation pressure 
acting on the gas counterbalances the gravitational at- 
traction of the BH: 



LEdd = 1-^-10 



38 



M 



BH 



Mr. 



erg s 



(8) 



The obtained Mbh and the relative QSO Eddington 
ratios are listed in Tab. IH The two relations, Eq. |4]and 
Eq. O lead to a difference in the mass estimates of a 
factor ~ 1.7. The results obtained for the sample of 
published sources via Eq. [5] are in good agreement with 
previous estimates, indicating that different fitting pro- 
cedures imply variations smaller than the errors due to 
the intrinsic scatter of the estimator. The average dif- 
ference between the BH masses estimated by us and the 
ones previously published is equal to 0.1 dex. OnJ^^in 
the case of the QSO SDSS J0005-0006 ( Kurk et al.ll2007l 
online Appendix Fig. A. 11) this difference is higher (0.7 
dex): in this case the Mgll line is severely affected by 
the atmospheric absorption and the line fit becomes sig- 
nificantly procedure-dependent. 

4.2. Dependence on the S/N and on the systematics in 
the observations 

In Fig. 14.21 we show the Mbh obtained via Eq. |4] as a 
function of the S/N per A of the Mgll peak. The noise 
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TABLE 3 

Fitted spectral properties. The spectral decomposition has been performed 
on the specific fluxes fx in units of 10"^'^ erg s~^ cm~^ a~^. (1) qso name. 
For sources with multiple observations: *, IJiang et al.I ()2Q07I) : + JKurk et al.I 
(2007).; (2) POWER-LAW slope; (3) Feii normalization; (4) Men central 
wavelength [a]; (5) fwhm [km s"-*^]; (6) measured specific flux of the 

continuum at 3000 a, in units of 10"^'^ erg s~^ cm~^ a~^. j0303 - 0019: the 
Feii major features are severely affected by the telluric absorptions, 

RESULTING IN A WEAKER CONSTRAINT ON THE FE II TEMPLATE NORMALIZATION WHICH IS 
IN ANY CASE CONSISTENT WITH 0. 



QSO name 


slope 


Fe II norm 


Me: II X^^^u 

xvig, 11 '""peak, 


Mgii FWHM 


Fx (3000A) 


(1) 




(2) 


(3) 








BR 1033 - 0327 


-2.4 


0.84+0.02 


2804.6+2.3 


4248+120 


8.20+1.59 


BR 0019 - 1522 


-2.3 


1.01+0.03 


2802.6+0.6 


4317+120 


7.78+1.50 


BR 2237 - 0607 


-1.3 


1.57+0.36 


2796.3+0.6 


3811+163 


19.98+0.12 


SDSS J0310 


- 0014 


-1.5 


0.24+0.05 


2797.3+1.1 


4087+260 


2.83+0.05 


SDSS J1021 


- 0309 


-1.7 


0.26+0.09 


2797.4+0.6 


3100+162 


2.76+0.09 


SDSS J0210 


- 0018 


-2.7 


0.11+0.06 


2798.4+2.1 


6543+803 


3.42+0.05 


SDSS J0211 


- 0009 


-2.4 


0.38+0.05 


2795.7+1.5 


5975+468 


1.44+0.28 


PC 1247 + 3406 


-2.2 


0.30+0.10 


2797.7+0.7 


4094+197 


3.86+0.08 


SDSS J0338 


+ 0021 


-1.8 


0.51+0.05 


2798.4+0.4 


2969+128 


3.20+0.05 


SDSS J1204 


- 0021 


-2.3 


0.83+0.03 


2796.2+3.0 


5753+208 


4.36+0.04 


SDSS J0005 


- 0006 


0.0 


0.18+0.04 


2796.2+0.1 


1036+65 


0.94+0.02 


SDSS J1411 


+ 1217* 


-1.7 


0.20+0.07 


2800.4+1.2 


2208+317 


2.39+0.06 


SDSS J1411 


+ 1217+ 


1.2 


0.36+0.04 


2787.8+0.4 


2824+168 


3.03+0.05 


SDSS J1306 


+ 0356* 


-2.3 


0.21+0.12 


2809.5+0.5 


3158+145 


1.92+0.05 


SDSS J1306 


+ 0356+ 


-2.5 


0.13+0.13 


2809.8+1.1 


4134+351 


1.97+0.04 


SDSS J1630 


+ 4012 


-0.7 


0.03+0.03 


2798.0+2.0 


3366+533 


1.46+0.30 


SDSS J0303 


- 0019 
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0.00+0.15 


2802.28+0.4 


2307+63 


0.60+0.01 


SDSS J1623 
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0.31+0.02 


2795.0+0.3 


3587+118 


1.81+0.05 


SDSS J1048 


+ 4637 


-1.7 


0.56+0.03 


2797.9+1.4 


3366+532 


4.22+0.04 


SDSS J1030 


+ 0524* 


-1.9 


0.16+0.07 


2807.2+0.6 


3449+151 


1.76+0.05 


SDSS J1030 


+ 0524+ 


-1.3 


0.17+0.08 


2805.9+0.9 


3704+231 


1.88+0.03 


SDSS J1148 


+ 5251 


-2.1 


0.91+0.01 


2793.7+0.7 


5352+289 


4.22+0.08 


SDSS J0353 


+ 0104 


-2.2 


0.43+0.06 


2799.6+0.7 


3682+281 


2.16+0.10 


SDSS J0842 


+ 1218 


-0.9 


0.38+0.06 


2794.2+0.7 


3931+257 


2.38+0.12 



was derived from regions of the spectra free of line emis- 
sion. There is no evidence of any dependence of the Mbh 
estimates on the S/N of the Mgll hne peak. For SDSS 
J1411 + 1217, SDSS J1306 + 0356 and SDSS J1030 + 0524 
we are in the fortunat e situation of havi ng two indepen- 
dent observation s bv [jiang et "aTl (l2QQ7l ). with Gemini- 
GNIRS, and by iKurk et al.I (j2007[ ). with VLT-ISAAC. 
We can thus use these to analyze the dependence of the 
Mbh estimates on the systematics in the observations. 
For these 3 sources Jiang et al. (2007) and Kurk et al. 
(|2007r ) derived (using Eq. [5]) 0.6, 1.1, 1.0 xlO^ and 
1.1, 2.4, 1.4 xlO^ M0, respectively. Using our fitting 
results and Eq. [5l we obtain 0.5, 0.9, 1.1 xlO^ M© and 
1.0, 1.7, 1.4 xlO^ respectively, in good agreement 
with the published values. For those objects whose red- 
shifts make their Mgll line fall on top of atmospheric 
absorption bands, the Mbh estimated from the two in- 
dependent spectra can be significantly different (up to 
0.3 dex in the case of SDSS J1306 + 0356): therefore 
we argue that the atmospheric contamination is a major 
contributor to systematics. 

4.3. Comparison with SDSS QSO sample at 
0.35 <z< 2.25 

To study the evolution of the BH population with cos- 
mic time, we compare our sources with a sam ple of SDSS 
QSO s at lower redshift (SDSS Data Release 7. lShen et al.I 
l2010f ). For objects with multiple observations we con- 
sider the weighted mean of the individual measurements. 
Since the Mbh estimates significantly depend on the 



adopted estimator (see lPetersQnll201Ql : iShen et al.ll2008f ). 
we select among the -100,000 QSOs in the SDSS DR7 
sample those with measurements of the Mgll FWHM 
and of ALa(3000A): i.e., a subset of '—47,000 sources 
with 0.35 < z < 2.25. We also consider a sample of nine 
z — 6 QSOs a t lower luminosities (M1450 < —24) from 
^Willo tt et al.I ([2010a'), with published measurements of 
the Mg II FWHM and of ALa(300 0A). We obt ain L^^/ us- 
ing the bolometric correction by iShen et aTl (j2008l) (see 
Eq. [7j) and we compute Mbh using Eq. HI in anal ogy 
with the analysis performed in our sample. In Fig. 14.31 
we plot the estimated Mbh as a function of the bolo- 
metric luminosity for the three samples. The black solid 
lines indicate regions of the parameter space with con- 
stant Eddington ratios: Li^oi/LEdd = 0.01, 0.1, 1.0. 
The 0.35 < z < 2.25 and our high-z QSO populations 
occupy different regions of the parameter space: the 
0.35 < z < 2.25 redshift population shows an average 
{log{Lboi/LEdd)) ^ -1.2 dex {Lboi/LEdd ^OM) with a 
dispersion of 0.3 dex, while the typical ratio for the high-z 
QSOs is {log{Lboi/LEdd)) ^ -0.3 dex {Lboi/LEdd ^0.5) 
with a d ispersion of 0.3 dex. The lower luminosity z ^ 6 
sample (jWillott et al.I [2QTQah is characterized by lower 
Mbh but comparable Eddington ratios than the most 
luminous SDSS QSOs at the same redshift. We there- 
fore conclude that the average Eddington ratio for the 
z ^ 6 QSO sample studied to date is significantly higher 
than the typical SDSS QSOs at 0.35 < z < 2.25 . 

We have to specify that IShen et "aTl (|2010l ) fit the 
spectra considering the power-law non-stellar contin- 
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Fig. 1. Reduced spectra of the QSOs SDSS J0842+1218 (upper-left panel), SDSS J0353+0104 (upper-right panel) and SDSS J2054-0005 
(lower panel). The observed spectra are shown in black continuous line. The modeled components are: power-law continuum (blue dotted 
line), Balmer pseudo continuum (purple dashed line). Fell normalized template (light blue dotted line), Mgll emission line (red dotted line). 
The sum of the first set of components (power-law continuum + Balmer pseudo continuum + Fe ll normalized template) is overplotted to 
the spectrum as a green solid line, while the sum of all the components is over plotted as a red solid line. The spectral decomposition was 
not performed for the weak line QSO SDSS J2054-0005. Telluric absorption bands are indicated over the spectra with the symbol 0: they 
are extracted from the ESO sky absorption spectrum measured on the Paranal site at a nominal airmass of 1. 
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Fig. 2. — The Fe ll normalization, obtained from the fit of the first 
set of components, depends on the power-law slope and its nor- 
malization (intercept). SDSS J0842+1218 (top panel) and SDSS 
J0353+0104 (bottom panel): domain analysis for Fell error 
computation. In each panel we show: a) two dimensional pro- 
jections of the 3D x^-surfaces (Fell normalization vs Intercept, 
upper-left plot; Fell normalization vs Slope, bottom-left plot; In- 
tercept vs slope, bottom-right plot): contours represent iso-x^ lev- 
els spaced by a factor of 2 while the best fit case is marked with 
a dot; b) probability distribution for the Fell template normal- 
ization (upper-right plot): the distribution has been obtained by 
marginalizing the 3-D probability distribution considering only the 
triplets for which x^ — Xmin < dashed vertical lines mark 

our estimate of the 1-cr confidence level. 

uu m, the Fell line forest (mod eled using the template 
by 'Vestergaard fc WilkesI "2001) , and both the narrow 
(simple Gaussian) and broad (simple/double Gaussian) 
components of the Mgll doublet. Since in our sam- 
ple we are not subtracting the narrow component of 
the Mgll doublet (not resolved), the mass estimates of 
the high-redshift QSOs may be systematically biased to- 
wards slightly lower values with respect to the ones in the 



0.35 < z < 2.25 redshift sample. Nevertheless, we note 
that our fitting proced ure yields r esults consistent with 
the ones published by J iang et "ahl ()2007[ ) who, except for 
the cases in which the line profile is severely affected by 
the atmospheric absorption, fits the Mgll emission line 
with a double Gaussian instead of a single one. There- 
fore, we conclude that the Mbh estimates do not strongly 
depend on the adopted fitting procedure. 



4.4. Comparison with luminosity matched SDSS QSO 

sample 

In order to perform a statistically robust comparison 
we now extract only objects with the same (high) bolo- 
metric luminosity: 10^^ < L})oi [erg s~^] < 3 x 10"^^ 
from our sample and from the SDSS comparison sam- 
ple at 0.35 < z < 2.25. In Fig. 1131 we plot the 
Mbh and the Edddington ratio as a function of red- 
shift. The average Mbh at high-z is systematically lower 
than the typical Mbh of the 0.35 < z < 2.25 popu- 
lation at similar luminosities, while the Eddington ra- 
tio is higher. Comparing the Eddington ratio distribu- 
tion of the two luminosity matched samples (Fig. 14. 4| 
10"^^ < Lboi [erg 5"^] < 3 x 10"^^) we obtain that for the 
0.35 < z < 2.25 population {log{Lboi/LEdd)) ^ -0.80 
dex {Li)oi/ LecLcL ~0.16) with a dispersion of 0.24 dex, 
while for the high-z one {log{Li)oi/ LEdd)) ^ —0.37 dex 
{Lijoi/LEdd ~0.43) with a dispersion of 0.20 dex. If we 
perform a Kolmogorov-Smirnov test on the two distribu- 
tions we obtain that the probability that they come from 
the same parent distribution is 7 x 10~^, i.e. negligible. 
This implies that the high-redshift QSOs in this luminos- 
ity bin are building up their mass more rapidly than the 
ones at z ^ 2 (in agreement with wh at has been found 
for subsam ples in previous studies, e.g. iKurk et ani2009l : 
iWillott et al. 2010). 

Another bias that could possibly affect our comparison 
between the high-z Eddington ratio distribution and the 
lower-z one is the intrinsic luminosity of the QSO broad 
lines. For instance the equivalent width of high ionization 
lines in AGN is kn own to anti-correlate with the con- 
tinuum luminosity (jBaldwinl[l977l : iBaskin fc Laorl 120041 : 



'■'plibli'sheidj z'<5.7. 
^Published, z~6 . 
• New data 



5 10 15 20 

S/N (per A) of the Mgll Peak 



Fig. 3. — Mbh estimates as a function of the S/N per A 
at t he Ms II peak. Blue squares: hterature sample, at z < 
5.7 (jlw amuro et al. 2002); black tr iangles: literature s ample , 
z ^ 6 (E arth et al. 2003; Iwamuro e t al. 2004; Jiang et all [20071 : 
IKurk et al.ii2007i . i200a ): red circles: this study. Our Mbh esti- 
mates do not depend on the S/N of the QSO spectra. 
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Fig. 4. — M-qh as a function of the bolometric luminosity: comparison between the present dataset and the SDSS QSOs sample at 
0.35 < z < 2.25. Bl ue square s: literat ure sample , z < 5 .7 (Iwamuro et al. 2002); black triangles: literature sample, z ~ 6 (Barth e t al.l 
[20031 : [Twamuro et al .' 2004; Jians_et^^2007; Kurk et al.l [2007. 2009); red diamonds: new ob servations; green pentagons: CFHQS QSO 
sample ( Willott et al. 2010); grey dots: SDSS QSOs sample 0.35 < z < 2.25 fDR7. IShen et al.ll2TC0l ): purple hues: isodensity contours. 
The black solid lines indicate regions of the parameter space with constant Eddington ratios: Li^oi/LEdd = 0.01, 0.1, 1.0. The 0.35 < 
z < 2.25 and our high-z QSO populations occupy different regions of the parameter space: the 0.35 < z < 2.25 redshift population shows 
{log{LijQi / L Edd)) ^ —1-2 dex with a dispersion of 0.3 dex, while the typical ratio for the high-z QSOs is {log^Lij^i / L ^dd)) ~ —0.3 dex with 
a dispersion of 0.3 dex. 

(SMBHs) in the early universe. The three outstand- 
ing theories for SMBH seed formation differ both in the 
astrophysical processes considered and in the resulting 
masses of the seeds: (i) BHs seeds of few hundreds of 
Mq can be produced by the first generation of stars 
(Pop III stars), formed out of zero metallicity gas; (ii) 
BH seeds of 10^ — 10^ Mq can be obtained via stellar- 
dynamical processes amongst Pop II stars; (iii) BH seeds 
of 10^ — 10^ can be produced via gas-dyna mical pro- 
cesse s (direct collapse of dense gas clouds) (see IVolonteril 
|2Q1Q| . for a review on the SMBH formation). 

We here consider only the 20 QSOs at z ~ 6, since 
they provide tighter constraints on the SMBH formation: 
11 sources f rom our sample {z > 5.7) and 9 QSOs by 
I Willott et al.l (L2010a). The resulting sample includes all 
the faint 2; ^ 6 QSOs known to date with estimated Mbh 
from the Mgll line. For this z ^ 6 population we obtain 
{Ioq^Mbh/Mq)) - 9.06 dex {Mbh ^ 1.1 x 10^ M©) with 
4.5. Impact on BH formation history a dispersion of 0.5 dex and {log{Lboi/LEdd)) ^ -0.23 

. ^ ^ r-n,TT 1 1 1 dex (L^o^/L^;^^^ '^0.6) with a dispersion of 0.29 dex. The 

Accurate measurements 01 BH masses and relative bd- ^ a a u -dxj a ^ ^ ^ ^ ^ 

^. ^ r^nr^ 1 • time uecdcd by a BH seed to grow at a constant rate 

dmgton ratios 01 z ~ QSOs can be used to give con- r • 1 i\/r ^ ^ ^ i\/r • 1 ^ 

^ ^. ^ ^1 r r • -n»TT irom au initial mass Mq to a final mass Mt is equal to 

stramts on the formation processes of super-massive BHs 



iShang et al.l[M)l : [Xueni '^20^ This effect is a func- 
tion of the line ionization potential and not observed in 
low ionization lines (e.g. Mgll line, Espey & Andreadis 
IT999). High-z QSOs are usually selected in function of 
their photometric colors as drop-out objects. This cri- 
terion is to some extent insensitive to the line bright- 
ness. For the SDSS QSOs sample with 0.35 < z < 2.25 
the existence of a broad line feature is instead required 
in addition to the color selection, resulting into a bias 
against weak line QSOs. Deca rli et al.l (120111) found that 
for the DR7 SDSS QSO sample (|Shen et al.l |2Q10) the 
Mgll line luminosity positively correlates with the con- 
tinuum luminosity. Selecting objects within a given Mg ll 
line luminosity bin would then result in a continuum lu- 
minosity cut, like the one previously taken into account. 
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Fig. 5. — Eddington ratio (upper panel) and Mbh (lower panel) 
as a function of redshift for subsamples of QSOs with bolomet- 
ric luminosities 10^^ < L^o^ [erg s"-*^] < 3 x 10^^ extracted 
from the 0.35 < z < 2.25 sample an d from our high-z sam- 
ple. Blue circles: z < 5.7 (llwamuro e t al. 2002); black circles: 
z 6 (Barth et ar'2003': Iwamuro et al. 2004: Jiang et al. 2007: 



IKurk et al. 2007, 2009, this study); grey circles: SDSS QSOs sam- 
ple 0.35 < z < 2.25 (DR7, Shen et al. 2010); purple hues: isoden- 
sity contours. The average Mbh at high-z is lower than the typical 
Mbh of the 0.35 < z < 2.25 population at similar luminosities, 
while the Eddington ratio is higher. 



0.6 




-1 -0.5 

log (kol/Ledd) 

Fig. 6. — Eddington ratio distributions for sub samples of QSOs 
with bolometric luminosities 10^^ < Ljjoi [erg s"-*^] < 3 x 10^^ 
extracted from the 0.35 < z < 2.25 sample and from our high-z 
sample: a) black dashed hue: SDSS QSOs at 0.3 < z < 2.2 (SDSS 
DR7, Shen et al. 2010); b) red solid hue: our sample at 4.7 < 
z < 6.5. The two distributions are intrinsically different: for the 
0.35 < z < 2.25 population {log{L}joi / L Edd)) ~ —0.80 dex with a 
dispersion of 0.24 dex, while for the high-z one {log(Liyoi / L Edd)) ~ 
—0.37 dex with a dispersion of 0.20 dex. The Kolmogorov-Smirnov 
probability that the two distributions come from the same parent 
distribution is negligible. 



(jVolonteri fc ReesI 120051 ): 



t = 0.45 Gyr x 



Ledd 1 f^t \ 



Lbol 



VMo 



where 0.45 Gyr is the Salpeter time, e is t he radiative effi- 
ciency (e ~ 0.07. IVolonteri fc Re esl l2005[ ). and Ltoi/Ledd 
is the Eddington ratio. A BH seed with mass Mq = 
10^, 10^ Mq accreting constantly at the Eddington ra- 
tio characteristic of our population {L}yoi/Ledd ^ 0-6), 
would then need a time t ~ 0.91, 0.53 Gyr, respectively, 
to obtain a final mass Mt equal to the mean Mbh of our 
sample. 

If the characteristic Eddington ratio of the ~ 6 QSO 
population is well represented by our mean estimate, 
models with BH seed masses ~ 10^ Mq become prob- 
lematic since the time needed for the seed to grow is 
comparable to the age of the universe at z ~ 6 (~ 0.9 
Gyr). Moreover, since our QSOs are selected from flux- 
limited surveys the resulting distribution of the Edding- 
ton ratios should be biased towards higher values: BHs 
with low accretion rates may n o t pass the QSQ selection 
magnitude limits. iKellv et al] ()2010l ) and IWillott et al.l 
(2010a) have shown that the intrinsic Eddington ratio 
distribution for a volume limited QSO sample is indeed 
shifted towards lower values with respect to a luminos- 
ity selected population. Our estimate hence represents 
an upper limit of the characteristic Eddington ratio of 
the 2: ~ 6 QSO population: lower Eddington ratios, ac- 
cording to Eq. m would yield even longer time for the 
BH seeds to grow. In the future, with complete sam- 
ples of QSOs at z ~ 6, it will be possible to give better 
constraints on the super massive BHs formation models. 

4.6. Fell/Mgll 

An estimate of the Fe/Mg abundance ratio at z ~ 
6 may serve as an indication of the onset of star 
formation in the highest-z QSOs (see Sec. [T]). The 
proxy usually used to trace the Fe/Mg abundance ra- 
tio at high-z is the Fell/Mgll line ratio. In the 
past, numerous NIR-spectroscopy studies have been car- 
ried out to analyze the Fe ll/M^ll line ratio in high-z 
QSO s (e.g. iMajoHno et al.li 2001. 2003; Penteric cTetall 
20021: llwamuro et al.' "2002, 2004; Barth et al] 12003 



Dietrich et al.l [2 003: Freudling et al. 2003; Willott et al.l 
20031 : iJiang etal.ii2007i : iKurk et al.h 2007). The combined 
results revealed an increase in the scatter of the mea- 
sured Fell/Mgll line ratios as a function of redshift (see 
iKurk et all l2007l ). A proposed explanation for the in- 
creased scatter is that some young objects have been ob- 
served such a short time after the initial starburst that 
the BLR has not been fully enriched with Fe yet. 

We compute the Fell flux by integrating the normal- 
ized Fe II template over the rest-frame wavelength range 
2200 A < Arest < 3090 A. The Mgii flux is computed 
by integrating the fitted single Gaussian over the range 
Xpeak i 5(7. The resulting Fell/Mgll flux ratios and rel- 
ative errors are listed in Tab. [H In Fig. 14.61 we show 
the evolution of the Fell/Mgll as a function of z (for ob- 
jects with multiple observations we consider the weighted 
mean of the indi vidual measureme nts). The reported 
scatter at 2; ~ 6 (|Kurk et al.ll2007l ) is now significantly 
reduced: our measurements of the line ratios span a 
range <Fell/Mgll< 6, while previous results were dis- 
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TABLE 4 

Estimated Mbh, QSO Eddington ratios and Feii/Mgii line ratios. The mass 
estimates have an accuracy of 0.4 dex that dominates the measurement 
uncertainties. (1) QSO name. For sources with multiple observations: *, 

IJlANG ET AL.I (|2Q07I ): + JKURK ET AL.I (|20Q7I ): (2) REDSHIFT ESTIMATE FROM MC II 
LINE; (3) Mbh IN UNITS OF 10^ Mq ESTIMATED USING EQ. (4] (4) Mbh IN UNITS OF 
10^ M0 ESTIMATED USING EQ. \E\ (5) EDDINGTON ACCRETION RATIO DERIVED FROM 

Mbh col. [2]; (6) Eddington accretion ratio derived from Mbh col. [3]; (7) 

Feii/Mgii ratio. 



QSO name 


z 


Mbh [10^ Mq] 


Lboi 


/ Ledd 


Feii/Mgii 








Eq.Hl 


Eq.[5r 


Eq.ll 


Eq.[5] 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


BR 1033 - 0327 


4.521+0.005 




Z.O 


0.5 


0.8 


2.6+0.1 


BR 0019 - 1522 


4.534+0.001 


A A 
4.4 


9 Q 
Z.O 


0.5 


0.7 


2.7+0.1 


BR 2237 - 0607 


4.561+0.001 


0.4 


A n 

4.U 


1.0 


1.3 


2.1+0.5 


SDSS J0310 


-0014 


4.672+0.002 


Z.O 


1 A 
1.4 


0.3 


0.6 


2.7+0.6 


SDSS J1021 


- 0309 


4.703+0.001 


1.5 


0.8 


0.5 


1.0 


3.3+1.1 


SDSS J0210 


-0018 


4.715+0.004 


7.2 


4.1 


0.1 


0.2 


1.5+0.9 


SDSS J0211 


- 0009 


4.894+0.003 


4.1 


2.1 


0.1 


0.2 


4.2+0.6 


PC 1247 + 3406 


4.929+0.001 


3.1 


1.9 


0.4 


0.7 


2.1+0.7 


SDSS J0338 


+ 0021 


5.027+0.001 


1.5 


0.9 


0.7 


1.2 


5.0+0.6 


SDSS J1204 


- 0021 


5.090+0.007 


6.8 


4.1 


0.2 


0.4 


3.3+0.1 


SDSS J0005 


- 0006 


5.844+0.001 


0.1 


0.06 


3.6 


7.1 


4.7+1.1 


SDSS J1411 


+ 1217* 


5.854+0.003 


0.9 


0.5 


1.3 


2.2 


2.0+0.7 


SDSS J1411 


+ 1217+ 


5.903+0.001 


1.6 


1.0 


0.9 


1.5 


2.9+0.5 


SDSS J1306 


+ 0356* 


6.017+0.001 


1.7 


0.9 


0.6 


1.0 


4.8+2.7 


SDSS J1306 


+ 0356+ 


6.018+0.003 


2.9 


1.7 


0.3 


0.6 


-•-•^-1.5 


SDSS J1630 


+ 4012 


6.058+0.005 


1.7 


0.9 


0.5 


0.8 


0.5+0.5 


SDSS J0303 


-0019 


6.079+0.001 


0.5 


0.3 


0.6 


1.3 


0.0+5.8 


SDSS J1623 


+ 3112 


6.211+0.001 


2.2 


1.2 


0.5 


0.8 


4.3+3.3 


SDSS J1048 


+ 4637 


6.198+0.004 


6.0 


3.9 


0.4 


0.6 


4.1+0.3 


SDSS J1030 


+ 0524* 


6.302+0.001 


2.0 


1.1 


0.5 


0.9 


2.9+1.3 


SDSS J1030 


+ 0524+ 


6.299+0.002 


2.4 


1.4 


0.5 


0.8 


2.3+1.2 


SDSS J1148 


+ 5251 


6.407+0.002 


7.4 


4.9 


0.3 


0.5 


4.6+0.2 


SDSS J0353 


+ 0104 


6.072+0.002 


2.4 


1.4 


0.5 


0.8 


4.7+0.7 


SDSS J0842 


+ 1218 


6.069+0.002 


2.9 


1.7 


0.4 


0.7 


3.1+0.5 



tributed up to Fell/Mgiw 15. The Fell/Mgll line mea- 
surements do not show any clear dependence on the esti- 
mated S/N per A of the spectral continuum (see Fig. l4.6p . 
On the other hand, our flux ratio es timates are well in 
agreement with previous results by iKurk et all ()2QQ7l ) 
who performed an analysis similar to ours. We can con- 
clude that this measurement is significantly dependent 
on the performed analysis, rather than on the S/N of 
the spectra. This result immediately implies that we do 
not need to invoke different chemical evolutionary stages 
of the BLR gas. We do not see evidence for evolution 
of the estimated Fell/Mgll line ratio as a function of 
cosmic age for 4 < z < 6.5. Moreover, if we consider 
the SDSS QSO sample at 0.35 < z < 2.25 CShen et al.l 
[2010. see Sec. 14. 3p . and we compute the Fell/Mgll line 
ratio starting from the published Fe ll and Mg ll fiux mea- 
surements, we see no evolution of the line ratio also for 
0.35 < z < 2.25. However it is not possible to per- 
form a direct comparison between the two samples, since 
these line ratio measurements are highly dependent on 
the adopted fitting procedure. 

To constrain the onset of the star formation in these 
early QSOs, we would need to derive an estimate of the 
Fe/Mg abundance ratio from the measurements of the 
Fell/Mgll line ratio. Unfortunately, an accurate conver- 
sion cannot be performed since the detailed formation of 
the Fe ll b ump in QSO spect ra has not been fully under- 
stood yet. IWills et aTl ()1985l ) discussed the suitability of 
the Fell/Mgll line ratio as a proxy of the Fe/Mg abun- 
dance ratio, as the regions in which the two ions are 



produced and the radiative transfer of the two lines are 
both very different. Computing the Fell and Mgll line 
strength for a gas with cosmic abundances and using a 
wide range of hydrogen densities and ionization parame- 
ters, they predicted typical Fe ll/Mg ll line ratios between 
1.5 and 4.0 (this is indicated in Fig. 14.61 as a grey-shaded 
area), and attributed higher Fe ll/Mg ll values to an over- 
abundance of Fe with respect to Mg. However these 
computations were bas ed on a limited 70-level model of 
the Fe atom. Baldw in et al.l (|2004) used a 371-level Fe+ 
model to reproduce the observed Fell emission proper- 
ties in AGNs considering a large set of models of broad 
emission line region clouds. From their results, all the 
observed Fe II features can be reproduced only if: (i) the 
BLR is characterized by the presence of significant micro- 
turbolence; (ii) the Fe ll emitting gas has different prop- 
erties (density and/or temperature) with respect to the 
gas emitting other broad lines. In any case the strength 
of the Fe ll emission relative to the emission line of other 
ions (e.g the Fell/Mgll line ratio) depends as much on 
the Fe abundance as it does on other physical param- 
eters of the BLR (e.g. turbulence velocity), making it 
difficult to convert the observed line ratios in abundance 
ratios. Nevertheless, the study of the Fell/Mgll line ra- 
tio as a function of look-back time can be used to give 
constraints on the BLR chemical enrichment history, un- 
der the assumption that physical conditions of the BLR 
that determine the Fe ll emission are not evolving. 

The observed lack of evolution in the measured 
Fell/Mgll line ratio can then be explained with an early 
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chemical enrichment of the QSO host. I.e. the QSOs 
in our sample must have undergone a major episode of 
Fe enrichment in a few hundreds Myr before the cos- 
mic age at which th e y ha ve been observed (~ 0.9 Gyr). 
[Matteucci & Recchi' (2003) showed that for a massive el- 
liptical galaxy characterized by a very intense but short 
star-formation history, the typical timescale for the max- 
imum SN la rate can be as short as 0.3 Gyr. This im- 
plies that if the Fe in high-z QSO hosts is mainly pro- 
duced via SN la explosions, it would be possible to ob- 
serve fully enriched BLR at 2; ^ 6. On the other hand, 
iVenkatesan e t al.l (|200 4) pointed out that SNe la are not 
necessarily the main contributors to the Fe enrichment, 
and that stars with a present-day initial mass function 
are sufficient to produce the observed Fell/Mgll line ra- 
tios at z ~ 6. Fe could also be generated by Pop III 
stars: these very metal poor stars with typical masses 
M > 100 M0 might b e able to produce large amounts of 
Fe within a few Myr (Heger Wooslevll2002"[ ). 
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Fig. 7. — Evolution of the Feii/Mgii line ratio for z > 4. Black 
points: literature sample, z > 4.4; red points: new targets. Grey- 
shaded area: expected Fell/Mgll line ratio obtained by model- 
ing the Fe 11 and Mg ll line strength considering gas with cosmic 
abundance and a range of ionizing continua (Wills et al. 1985). 
The reported scatter at 2; ~ 6 (K urk et al.h2007i ) is now signifi- 
cantly reduced: our measurements of the line ratios span a range 
<Fell/Mgll< 6, while previous results were distributed up to 
Fell/Mgll~ 15. No redshift evolution of the Fell/Mgll ratio is 
observed for 4 < z < 6.5. 



5. SUMMARY AND CONCLUSIONS 

We presented NIR-spectroscopic observations of three 
z - 6 SDSS QSOs. Our NIR spectra cover the Mgll 
and Fell emission features, which are powerful probes 
of Mbh and of the chemical enrichment of the BLR. 
The new data extend the existing SDSS sample towards 
the faint end of the QSO luminosity function. We have 
collected 22 literature spectra (19 different sources) of 
high-redshift {z > 4) QSOs covering the rest-frame wave- 
length range 2700 A < A < 3200 A. The final sam- 
ple is composed of 22 sources: our three new sources 



_ "PuiDlislied' sample,' z<h.^ 
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6 "•New data 
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S/N (per A) of the continuum at 3000 A 

Fig. 8. — Feii/Mgii line ratio estimates as a function of the S/N 
per A of spectral continuum at 3000 A. Blue squares: literature 
sample, at z < 5.7 (Iwamuro et al. 2002); black triangles: literature 
sample, z ^ 6 (.Ba rth et al. 2003; Iwa muro et a l. 2004; Jia ng et all 
[20071 : [KuFk et aT.! 2007, 2009); red circles: new observations. The 
Fell/Mgll line measurements do not show a dependence on the 
S/N per A of the spectral continuum. 

at z ~ 6, 10 s pectra from the literature of QSOs with 
4 < z < 5.7 ([Iwamuro et alJ [2QQ2I) a n d 9 spectra of 
QSOs with z ^ 6 (|Barth et al] 12003': 'Iw amuro et al.l 
[2003 : lJiang et al.ll2007l : iKurk et al.l l2007. 2009^" In order 
to see how the derivation of physical parameters depend 
on the fitting strategy employed, we performed a new 
consistent analysis of the sample and we gave an estimate 
of the Mbh, of the Li^oi/LEdd and of the Fell/Mgll line 
ratio. 

Our results and conclusions can be summarized as fol- 
lows: 

1. We have estimated the Mbh from the Mgll emis- 
sion line using empirical mass scaling relations. 
The QSOs in our sample host BHs with masses of 
^ 10^ M0. Our results are in good agreement with 
previous estimates, indicating that different fitting 
procedures imply variations smaller than the errors 
due to the intrinsic scatter of the Mbh estimator. 

2. High-redshift QSOs are accreting close to the Ed- 
dington luminosity: {log{Li)oi/ LEdd)) ^ —0.3 dex 
{Ljjoi/LEdd ^0.5) with a dispersion of 0.3 dex. The 
distribution of observed Eddington ratios is sig- 
nificantly different (average Eddington ratio ^ 8 
times higher) from that of a comparison sample 
at 0.35 < z < 2.25 (SDSS, DR7). This result is 
not biased by the luminosity selection of the high-z 
sample: the difference in the two Eddington ratio 
distributions persists even considering two luminos- 
ity matched subsamples {z ^ 6 average Eddington 
ratio ~ 3 higher than lower redshift one). The 
high-z sources are accreting faster than the ones at 
0.35 < z < 2.25. 

3. We calculated fluxes for the Mgll and Fell lines. 
The obtained values are not always in agreement 
with the published ones, indicating that these 
measurements are dependent on the performed 
analysis. The previously observed scatter in the 
Fell/Mgll line ratio at 2; ~ 6 is significantly re- 
duced in this work: the BLRs in the highest-z 
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QSOs studied to date all show comparable level 
of chemical emichment. 

4. No redshift evolution of the Fell/Mgll ratio is 
observed for 4 < 2: < 6.5. If we consider the 
Fell/Mgll line ratio as a second order proxy for 
Fe/Mg, this indicates that the QSOs in our sample 
must have undergone a major episode of Fe enrich- 
ment in the few hundreds Myr preceding the cosmic 
time at which the sources are observed. 

5. Our results do not show any significant dependence 
on the S/N of the spectra: the Mbh estimates do 
not correlate with the S/N per A of the Mgll line 
and the Fe ll/Mg ll line ratios do not show any clear 
dependence on the S/N per A of the spectral con- 
tinuum. 

Even though the currently faintest known QSOs are 
included in the study we presented here ( Willott et al . 



l2QlQal : [Kurk et al.ll2QQ9'l ). we are still not probing the bulk 
of QSOs that should be present at that redshift. Future 
studies of QSOs that populate the faint end of the QSO 
luminosity function at z ~ 6 are needed to investigate 
whether or not the results found here are applicable to 
ah ^ - 6 QSOs. 
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APPENDIX 
PIT RESULTS: LITERATURE DATA 

Hereafter we show the fit results for the Hterature sample and the relative maps for Fell error computation 
(analogous to Fig. 13. 4p . The sources are sorted in redshift. 

In the right panels we show the spectral decomposition. The observed spectra are shown as a black continous line. 
The modeled components are: power-law continuum (blue dotted line), Balmer pseudo continuum (purple dashed 
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line), Fell normalized template (light blue dotted line), Mgll emission line (red dotted line). The sum of the first set 
of components (power-law continuum + Balmer pseudo continuum + Fe ll normalized template) is overplotted to the 
spectrum as green solid line, while the sum of all the components is overplotted as a red solid line. Telluric absorption 
bands are indicated over the spectra with the symbol 0: they are extracted from the ESO sky absorption spectrum 
measured on the Paranal site at a nominal airmass of 1. 

The Fell normalization, obtained from the fit of the first set of components, depends on the power-law slope and 
its normalization (intercept). In the left panel we show the domain analysis for Fell error computation: a) two 
dimensional projections of the 3D x^-surfaces (Fell normalization vs Intercept, upper-left plot; Fell normalization vs 
Slope, bottom- left plot; Intercept vs slope, bottom-right plot): contours represent iso-x^ levels spaced by a factor of 2 
while the best fit case is marked with a dot; b) probability distribution for the Fell template normalization (upper- 
right plot): the distribution has been obtained by marginalizing the 3-D probability distribution considering only the 
triplets for which x^ — Xmin ^ 1' dashed vertical lines mark our estimate of the 1 — a confidence level. 
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